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Pui EFFICIENT METHOD FOR 


CARTESIAN TRANSFORMATIONS 
By Saiiuel A. Pines and Paul F. Flanagan 


1.0 SUMMARY AND INTRODUCTION 


In this note, an effirient alternate method used to perform 
Cartesian coordinate transformations is derived. The nethod replaces 
the 3 by 3 matrix product transformation procedure by a simpler vector 
operation which requires less computer storeme and less computer 
execution time. Also, for applications in which interpolation routines 
€ 0 *e applied to stored transformation matrices , the method outlined in this 
document is more efficient because it requires the interpolation and storage 
of only four elements, which results in a more accurate final computation. 


2.0 CARTESIAN COORDINATE TRANSFORMATIONS BY MEANS 
OF RIGID VECTOR ROTATIONS 


A 3 by 3 transformation matrix A is given which causes the 
Cartesian coordinate system I to be transformed to the Cartesian 
coordinate system II. 

Thus, any vector R^ in the original, coordinate system is carried 
into the second coordinate system by 

Rjj * ARj (1) 

The same transformation can be achieved by a rigid rotation C'f 
Rj about a unit vector N through some euxgle w [eq. (2)]. The 

resultant vector equation may be found in reference 1. 

Rjj = Rj cos w N»Rj( 1 - cos w)N + sin w NxR^ 


( 2 ) 


2 


If the matrix is required, it can be generated from the unit vector 
N and the angle w by the following equation. 


A 


I cob w (l - cos w)NN 


T 


sin w Nx 


(3) 


where 


Nx « 


0 






0 


The vector N and the trigonometric function cos w are cosputed 
from the matrix as follows. 


cos w 


trace A - 1 


sin w “ cos^w 


(Ua) 


(Ub) 


N » [(a32 - ag3)^ + (a^_3 - a3^)^ + - a^g)^] 


- 1/2 


32 

'l3 


‘21 


- a 


- a 


23 

31 


- a 


12 


ihc) 


3.0 DERIVATION OF THE EQUIVALENT TRANSFORMATION EQUATION 


In reference 1, it is shown that the rigid rotation matrix i.*.ay 
be written as the sum of a pymmetric and a skew symmetric matrix. 


The symmetric matrix terms of equation (3) are 


( 5a) 


A sym * I cos w + (l - cos w) NN 
The skew symmetric terms are given by 

A skew sym = sin w Nx (5b) 

The trace of A is given by the trace of its symmetric part; thus, 
trace A - 3 cos w + (l - cos w)(N^^ + N^^ + N^^) 


cos w 


' 6 ) 


3 


Because N 


is a unit vector, equation (7) 


is obtedned. 


cos w ■ 


trace A - 
2 


1 


( 7 ) 


Because a rigid rotation about N through an angle w greater 
than 7T may be replaced by a rigid negatiye rotation about the name 

vector N through 2tt - w, the angle w can be restricted to the 
first or second quacrants, and 

sin w ■ ^1 - cos^ w (8) 


To obtain the components of the vector N , note that the sk »w 
symmetric part of A is given by 


A skew sym 



(9) 


This expression must be equal to the skew symmetric paurt of 
A (N, w) so that 



a'^ 

— » sin V Nx 


( 10 ) 


It follows that because N is a unit vector, equation U(c) 
is obtained. 


I 


I 
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